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We study the dynamics of a ferromagnetic domain wall driven by an external magnetic field 
through a disordered medium. The avalanche- like motion of the domain walls between pinned 
configurations produces a noise known as the Barkhausen effect. We discuss experimental results on 
soft ferromagnetic materials, with reference to the domain structure and the sample geometry, and 
report Barkhausen noise measurements on Fe2iCo64Bi5 amorphous alloy. We construct an equation 
of motion for a fiexible domain wall, which displays a depinning transition as the field is increased. 
The long-range dipolar interactions are shown to set the upper critical dimension to dc = 3, which 
implies that mean-field exponents (with possible logarithmic correction) are expected to describe the 
Barkhausen effect. We introduce a mean-field infinite-range model and show that it is equivalent to 
a previously introduced single-degree-of-freedom model, known to reproduce several experimental 
results. We numerically simulate the equation in d = 3, confirming the theoretical predictions. 
We compute the avalanche distributions as a function of the field driving rate and the intensity 
of the demagnetizing field. The scaling exponents change linearly with the driving rate, while the 
cutoff of the distribution is determined by the demagnetizing field, in remarkable agreement with 
experiments. 

PACS numbers: 75.60.Ej, 75.60.Ch, 64.60.Lx, 68.35.Ct 



I. INTRODUCTION 

The Barkhausen effect vi^as first observed in 1919 
recording the noise produced by the sudden reversal 
of Weiss domains in a ferromagnet. Since then, the 
Barkhausen effect has been vsridely used as a non- 
destructive method to test magnetic materials and a de- 
tailed statistical analysis of the noise properties has been 
performed [|[D. Beside its practical and technological 
applications, the Barkhausen effect has recently attracted 
a growing interest as an example of a complex dynamical 
system displaying scaling behavior. It has been experi- 
mentally observed that a histogram of Barkhausen jump 
sizes follovirs a power law , a result which has analo- 
gies with other driven disordered systems, ranging from 
flux lines in type II superconductors to microfractures 
and earthquakes |l^ , where the dynamics takes place 
in avalanches. While the ambitious goal to build a com- 
mon theoretical framework for all these phenomena is 
still far from being reached, theoretical analysis of each 
system might shed light on the entire issue. 

In the case of the Barkhausen effect, the task is to 
explain the statistical properties of the noise, such as 
jump size distributions and power spectra, in terms of 
the microscopic details of the magnetization process. In 
general, three different mechanisms are involved during 
the process [pd| : domain nucleation and coalescence, co- 
herent spin rotation, and domain wall motion. Their 
different relevance along the hysteresis loop is in general 
very complicated and not easily predictable, as it depends 
on material properties, annealing conditions, and the ge- 
ometry of the sample. The Barkhausen noise is mainly 
due to the domain wall motion, therefore it is customary 



to study soft magnetic materials where a well defined 
domain structure is present and coherent spin rotation 
does not take place: in this case, once the structure is 
formed, the magnetization process takes place by motion 
of domain walls, rather than nucleation of new domains, 
which has a higher energetic cost due to magnetostatic 
interactions. 

The classical theoretical approach to the problem fo- 
cuses on the motion of the domain walls and their interac- 
tion with the disorder present in the medium. The simple 
schematization of the domain wall as a point moving in 
a random pinning field |l^ has been successfully used 
in the past to explain several properties of ferromagnetic 
materials, such as the Rayleigh law |l^. A theoretical 
analysis of the Barkhausen effect has been carried out 
in the same spirit iQ . Most of the measured properties 
can be reproduced by the model proposed by Alessan- 
dro, Beatrice Bertotti and Montorsi (ABBM) [|l5|. The 
crucial hypothesis of this model is that the pinning field 
is a random walk in space. This assumption is consistent 
with experiments but its microscopic justification is 
still unclear. In fact, an estimate ||l6| of the correlation 
length of the impurities typically present in the mate- 
rial gives a value much smaller than the one employed in 
jisf , implying that a Brownian pinning field can only be 
considered to be an effective picture. 

Recently, Urbach et. al. have proposed relating 

the properties of the Barkhausen effect to the depinning 
transition of an elastic surface in a random medium, a 
topic that has been studied extensively in recent years 
p9[ . The comparison between the values of the expo- 
nents predicted for the depinning transition and most 
experimental data was however unsatisfactory. 
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A completely different approach has been undertaken 
by Sethna et al. [^0|-^, who study field-driven nu- 
cleation in a non-equilibrium random-field Ising model 
(RFIM). In this model domain nucleation and growth 
are treated in the same way. When the external field 
is increased from the negative saturation, the spins flip 
to align with the local magnetization, eventually causing 
avalanches of neighboring spins. A thorough investiga- 
tion of this model shows that there is a second order 
critical point controlled by the amplitude of the disor- 
der . The power law distributions of the Barkhausen 
noise would then be related to the proximity of this crit- 
ical point [Q. This model neglects dipolar interactions 
and demagnetizing effects which are known to play a cru- 
cial role in the formations of domains, so its applicability 
to most experimental situations seems questionable. 

Here we approach the problem studying the motion of a 
flexible domain wall driven through a disordered medium. 
One of our aims is to bridge the gap between "classical" 
approaches to ferromagnetism ^2 13 1 and modern theo- 
ries of surface growth in disordered media ||l9|. In this 
way, we are able to clarify several assumptions present in 
phenomenological models of domain wall dynamics and 
to understand their limitations. 

We consider the case of an anisotropic material mag- 
netized along the easy axis, with 180° domain walls sep- 
arating regions of opposite magnetization (Fig. The 
disorder, due for example to non-magnetic inclusions or 
residual stresses, pins the domain wall motion which is 
driven by the external magnetic field. We assume that 
the disorder is localized and is either uncorrelated in 
space, or is only short-range correlated. The domain wall 
is assumed to be flexible, the stiffness beingdue to ferro- 
magnetic and magnetostatic interactions ||l^ , |2^ , p^ , and 
can therefore deform because of the local configurations 
of the disorder. The resulting equation of motion is dif- 
ferent from the one proposed by Urbach et al. ||l^], who 
treated incompletely dipolar interactions. Narayan fl^ ] 
has also considered dipolar interactions in this context, 
but his approximate analysis does not apply to d = 3 — 
the physical dimension for most of the experiments. 

We shall find that the scaling properties of the 
Barkhausen noise arise from the critical behavior ex- 
pected close to a depinning transition. The dipolar in- 
teractions generate a long-range term in the equation of 
motion which reduces the upper critical dimension from 
dc = 5, obtained for elastic interfaces [^,^, to dc = 3. 
Indeed, we shall see that mean-field critical exponents 
describe quite well a large amount of experimental data. 

The geometry of the sample has an important effect 
on the experimental results. A true depinning transition 
can only be observed when demagnetizing effects, oppos- 
ing the motion of the wall, are absent or very small. Oth- 
erwise, when the external field is increased at a constant 
rate, the wall is driven to a stationary motion around 
the depinning transition. The scaling is controlled by the 
external field driving rate and by the intensity of the de- 
magnetizing field, which in general depends on the shape 



of the sample. In particular, the driving field determines 
the exponents of the jump distributions while the cutoff 
is controlled by the demagnetizing field. 



We first introduce a mean-field interface model, in 
which the interaction range is infinite. Since the upper 
critical dimension is dc = 3, we expect that its criti- 
cal properties should agree with the three-dimensional 
model. Interestingly, we find the infinite-range model 
to be equivalent to the ABBM model. This observation 
explains why the ABBM model works so well in describ- 
ing the experimental data: it provides an effective one- 
degree-of-freedom description of the complex motion of a 
flexible interface. The elastic interactions along the wall 
moving in an uncorrelated medium give rise to an effective 
correlated pinning field experienced by the center of mass 
of the wall. In other words, the long-range correlations in 
the effective pinning field are not due to the correlation 
in the impurities present in the material. We note that a 
similar idea underlies the variational replica approach for 
equilibrium elastic interfaces in random media p^ ] , where 
one describes the complicated interactions between many 
degrees of freedom of the interface as a single particle in 
an effective potential. 



Finally, we simulate the full three-dimensional inter- 
face model and confirm the value of the upper critical 
dimension. We find that the results on three dimensional 
model do not fully agree with the mean-field predictions. 
In particular, the correct scaling of the cut-off can not be 
predicted by the infinite-range and ABBM models. The 
results of the simulations, however, agree remarkably well 
with experiments. 



The paper is organized as follows: In section II we 
discuss the experiments on the Barkhausen effect, intro- 
ducing the various scaling exponents. We briefly report 
experiments on an as cast Fe2iCo64Bi5 amorphous al- 
loy. In section III we construct the equation of motion 
for the dynamics of the domain wall. In section IV we 
derive the upper critical dimension and the mean-field 
exponents. In section V we derive scaling relations be- 
tween the critical exponents. In section VI we study the 
dynamics of the infinite range model as a function of 
the driving rate and the demagnetizing field. In section 
VII we present the result of numerical simulations. Sec- 
tion VIII is devoted to conclusions and discussion of open 
problems. A brief report of a subset of these results ap- 
pears in Ref. [p8[. 
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FIG. 1. The domain structure of Fe64Co2iBi5 amorphous 
alloy observed using the Kerr efltect. The domains are sep- 
arated by walls parallel to the magnetization. This is the 
typical structure observed in soft ferromagnetic materials. 



II. EXPERIMENTAL RESULTS 

The experimental results on the Barkhausen effect 
form an enormous body of literature that spans almost 
the entire century , but precise experimental results 
for the statistics of Barkhausen jumps have been reported 
only recently [^-Q . The distribution of Barkhausen jump 
sizes, measured at low driving rates, shows typically a 
power law behavior, but the scaling exponents reported 
in the literature span a wide range of values |^ . For this 
reason, it is important to carefully discuss the various 
experimental conditions, material properties and statis- 
tical uncertainties before direct comparison with a theory 
could be made. 

Under well-defined experimental conditions the results 
show a remarkable degree of universality: the scaling ex- 
ponents do not depend on the particular sample used 
P,p p5|j29| p^ . The measurements are taken only in the 
central part of the hysteresis loop around the coercive 
field, where domain wall motion is dominant while do- 
main nucleation and coherent spin rotations are negligi- 
ble 1 15| . The typical domain structure observed in these 
conditions is reported in Fig. |l|. Experiments were per- 
formed using a triangular waveform for the external field 
and different driving rates were employed. 

The signal amplitude distribution, directly related to 



the domain wall velocity, decays as power law |||,^,^,|o| 

P{v) - exp{-v/vo), (1) 

where c is proportional to the field driving rate and vq 
is the value of the cutoff. The avalanche size s (the area 
under the jump) and duration T distributions also decay 
as power laws and are very well fitted by (5l0,|29|] 



P{s) ^ s--f{s/so) T 
PiT) ^ T-''g{T/To) 



3/2- c/2. 



(2) 
(3) 



These laws have been tested for a variety of materials, 
such as amorphous (Co-base and Fe-base) ^9 33 and al- 
loys (Fe-) 1^,^. In Fig. ^, we report the size and duration 
distributions measured in an as cast Fe2iCo64Bi5 amor- 
phous alloy for different field driving rates. The experi- 
ments have been performed using the setup described in 
Ref. dj. The exponents agree perfectly with Eqs. (|[]|). 

The dependence of the exponents on the field driving 
rate can explain the variability in the experimen- 
tal values reported in literature, since many experiments 
were performed using a single linear driving rate or a 
sinusoidal one. Moreover, one should also be aware that 
the properties of the noise and thus the scaling exponents 
and the cutoff can change considerably through the hys- 
teresis loop ]l5| , ^ when domain nucleation and coherent 
spin rotations become relevant. 

To test the effect of the demagnetizing field, we per- 
form experiments on strips of with different lengths of an 
as cast Fe2iCo64Bi5 amorphous alloy. The intensity of 
the demagnetizing field decreases for longer samples. We 
find that the cutoff of distributions scales as sq ~ 1/fc 
and To ^ (Fig- H), where k is proportional to the 

intensity of the demagnetizing field (see Section III). We 
obtain the same results controlling k by changing the air 
gap between the sample and a magnetic joke. A com- 
plete account of these experiments will be deferred to a 
forthcoming publication. 

The power spectrum S(f) of the noise does not show 
in general such a marked robustness and is not described 
by a frequency- independent exponent: at low frequency 
/, 



S{f) ~ 



(4) 



where tp varies between ijj ~ 0.6 in Fe-Si, to ip I in 
amorphous alloys [ p9pl J3^ . After a crossover frequency, 
that depends on c, it decays with an exponent varying 
between -1.6 and -2 |6|j|,|2|,|l|,|| . When only a single 
domain wall is present the power spectrum was found to 
decay as /^^ Moreover, it has been observed that the 
power spectrum amplitude scales linearly with c. From 
the point of view of applications, it is important to distin- 
guish universal properties from material-dependent prop- 
erties that could be relevant to characterize the sample. 

In toroidal or frame geometries the demagnetizing field 
is practically absent and the magnetization process is 
quite different from in the previous case. The hystere- 
sis loop, instead of showing a extended linear part with 
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a stationary Barkhausen signal, displays a square form 
with a huge Barkhausen jump: the domain walls undergo 
a depinning transition as a function of the field. When 
the external field H exceeds the coercive field He, the do- 
main walls start to move with a velocity v that typically 
scales linearly with the field 
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FIG. 3. a) Distribution of Barkhausen jump durations 
measured in Fe2iCo64Bi5 amorphous alloy for different sam- 
ple lengths, b) The data collapse shows that To ~ 



This law has been first observed about 50 years ago by 
Williams, Shockley and Kittel |^ in a single crystal 
Fe-Si frame, and later confirmed for a variety of other 
soft ferromagnetic materials |^^. Before the onset of 
collective domain wall motion, one observes a series of 
Barkhausen jumps of increasing amplitude |^], but to our 
knowledge a quantitative analysis in terms of scaling ex- 
ponents has never been reported. 



FIG. 2. Distributions of Barkhausen jump sizes (a) and 
durations (b) measured in an as-cast Fe2iCo64Bi5 amorphous 
alloy for different driving rates. The lines are the fit with 
r = 3/2 — c/2 and a — 2 — c. The distributions have been 
obtained recording 6 ■ 10^ avalanches. 



III. DOMAIN WALL DYNAMICS 

The first thermodynamic theory of ferromagnetic do- 
mains is due to Landau and Lifshitz | |38[ |, who explained 
the presence of domains by energetic considerations. In a 
uniformly magnetized specimen, the discontinuity of the 
normal component of magnetization across the boundary 
of the sample creates a field that raises the total energy 
of the system. The creation of domains decreases this en- 
ergetic contribution at the price of a higher cost in wall 
energy. One can obtain a rough estimate of the number 
of domains by simply balancing these two terms. 

In order to describe accurately the magnetization pro- 
cess, it is necessary to analyze in detail the interactions 
present. In most ferromagnetic materials, due to the 
magnetocrystalline anisotropy or to the shape of the sam- 
ple, the magnetization has preferred directions. In the 
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simplest situation, there is a single easy axis of magneti- 
zation and the domains are separated by surfaces parallel 
to the magnetization, spanning the sample from end to 
end (see Fig. |l|). The domain walls are in general flex- 
ible, since local inhomogcneities can impose distortions 
of the surface, which would be flat in a perfectly ordered 
system. In some particular geometry in which demagne- 
tizing effects are minimized, it is even possible to obtain 
a single domain wall Q . 

We study the dynamics of a single 180° domain wall 
separating two regions with opposite saturation magne- 
tizations, directed along the x axis. If the surface has no 
overhangs, we can describe the position of the domain 
wall by a function h{f,t) of space and time (see Fig. ^). 
The equation of motion for the wall is given by 



^dh{r,t) _ dE{{h{r,t)}) 



dt 



Sh{f, t) 



(6) 



where E{{h{r,t)}) is the total energy functional for a 
given configuration of the surface and F is an effective 
viscosity. The motion of the domain wall is overdamped, 
since eddy currents cancel inertial effects, and thermal 
effects are negligible. 

We can split the energy into the sum of different con- 
tributions due to magnetostatic and dipolar fields, ferro- 
magnetic and magnetocrystalline interactions, and dis- 
order. In the following, we will express the energy in IS 
units. 



A. Magnetostatic fields 

In the presence of an external field H along the easy 
axis of magnetization the magnetostatic energy of the 
system is given by 



(7) 



where Ms is the saturation magnetization per unit vol- 
ume. 

Another contribution to the magnetostatic energy 
comes from the discontinuity of the normal component 
of the magnetization across the boundary of the sample. 
This generates an effective magnetic field, the so called 
demagnetizing field, that is opposed to the direction of 
the total magnetization. In some particular geometries 
(e.q. a uniformly magnetized ellipsoid) this field is con- 
stant along the sample. For a generic domain structure, 
an explicit expression for the demagnetizing field is often 
not available, but we expect in first approximation that 
the intensity of the demagnetizing field will be propor- 
tional to the total magnetization. Considering the field 
constant through the sample, its energy can be written 
as 



E, 



V 



(fr h{r,t) 



(8) 



where the demagnetizing factor Af takes into account the 
geometry of the domain structure and the shape of the 
sample and V is the sample volume. This term was also 
considered by Urbach et al. ||l^. The demagnetizing ef- 
fect can be avoided in suitable geometries, as in frame 
or toroidal specimens, but it is present in many common 
experimental situations. 



B. Dipolar interactions 

An effect similar to the one discussed above takes place 
inside the sample, where the local curvature of the surface 
can in general give rise to discontinuities in the normal 
component of the magnetization. We treat this effect 
introducing a "magnetic charge" density, which for a do- 
main wall separating two regions of magnetizations Mi 
and M2 is given by |39| 



cr = (Ml - M2) ■ fl 



(9) 



where n is normal to the surface. This charge is zero only 
when the magnetization is parallel to the wall. For small 
bending of the surface, we can express the charge as (see 
Fig. I) 



a(r) = 2MsCOs6l ~ 2M, 



dh{r, t) 
dx 



(10) 



where 9 is the local angle between the vector normal to 
the surface and the magnetization. The energy associ- 
ated with a distribution of charges a is given by 



OTT 



(11) 



Inserting the expression for a in Eq. ( |llD and integrating 
twice by parts, we obtain 



Ed^ / (Prcfr'h{f,t)K{f-f')h{r',t) 



(12) 



where the non local kernel has the form |23l 



27r|r - f'|3 y 1?^- f'p 



The interaction is long range and anisotropic, as can be 
seen by considering the Fourier transform 



Kip,q) 



(14) 



where p and q are the two components of the Fourier 
vector. 

In the preceding derivation we have implicitly assumed 
that the medium is infinitely anisotropic, so that the 
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magnetization never deviates from the easy axis. In prac- 
tice, however, the magnetization will rotate slightly from 
the easy axis because of the field created by the surface 
charges. A local change in the magnetization produces 
additional volume charges whose density is given by 



p(f) = V • M. 



(15) 



Neel has explicitly treated this effect obtaining an 
expression for the energy in the form of Eq. (|l2|) with a 
modified kernel 



effect of the disorder by introducing a random poten- 
tial V{f, h), whose derivative gives the local pinning field 
rj{f, h) acting on the surface. We take this random field 
to be Gaussian-distributed and short range-correlated 



(7/(f, h)T]{f", h')) = 5\r- r')R{h - h') 



(20) 



where R{x) decays very rapidly for large values of the 
argument. 



k{p,q) 



(16) 



where Q is a material dependent constant, whose value 
ranges from 5 to 10. This calculation shows that the qual- 
itative features of the interaction do not change if a finite 
magnetocrystalline anisotropy is taken into account. 

For the analysis we will perform later, it is important 
to generalize the kernel in any dimension. It is straight- 
forward to show that for d > 3 the kernel scales as 



K{q) cx 



(17) 



where and q± are the components of q parallel and 
perpendicular to the magnetization. 



C. Surface tension and disorder 

The magnetocrystalline and exchange interactions are 
responsible for the microscopic energy associated with 
the domain wall. While a very sharp change of the 
spin orientation has a high cost in exchange energy, a 
very smooth rotation of the spins between two domains 
is prevented by the magnetocrystalline anisotropy. The 
balance between these two contributions determines the 
width of the domain wall and its surface energy. The 
total energy due to these contributions is proportional to 
the area of the domain wall 



Edru = 1^0 I dVVl+ |V/l(f,i)|2 



(18) 



where vq is the surface tension. Expanding this term for 
small gradients we obtain 



Edw — 1^0 S, 



dw 



J d\\Wh{f,t)\\ 



(19) 



where Sdw is the domain wall area. This is the typical 
term associated with elastic interfaces. 

The disorder present in the material in the form of 
non-magnetic impurities, lattice dislocations or residual 
stresses is the reason for the jumps in the magnetization 
curve and for its hysteretic behavior. We can treat the 



D. The equation of motion 



Collecting all the energetic contributions, we obtain 
the equation of motion for the domain wall [ p8| . In order 
to avoid a cumbersome notation, wc will absorb all the 
unnecessary factors in the definitions of the parameters. 
The equation then becomes 



dh{r,t) 
dt 



= H ~kh + i/oV^/i(f, t) 



(frK{f- f'){h{r') - h[r)) + r^{f, h), 



(21) 



where the kernel K is given by Eq. (|l3|), k = A^oAfM^ 
ft, = / (Pr'h{r' , t)/V. Apart from the non-local ker- 
nel, this equation is similar to the equation proposed by 
Urbach et al. , which on its turn reduces when fc = to an 
elastic interface driven in quenched disorder. When the 
field is slowly increased, the demagnetizing field provides 
a restoring force that keeps the motion around the de- 
pinning transition. As we will show later, the non-local 
kernel changes the upper critical dimension, and hence 
the exponents, from the case of the clastic interface. 



h 



-M 




X 



FIG. 4. A domain wall separating two regions of opposite 
magnetization. The discontinuities of the normal component 
of the magnetization across the domain wall produce magnetic 
charges. 
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IV. MEAN-FIELD THEORY AND UPPER 
CRITICAL DIMENSION 

The mean-field theory provides a good quahtative de- 
scription of the depinning transition ]4l|-|4^. We will 
consider first the case fc = and H constant, which cor- 
respond to a conventional depinning transition. We will 
discuss in section VI the case k > 0, H (x t. Here, we 
proceed as in Refs. |2^:^^ , considering an infinite ranged 
interaction kernel in the equation of motion. To this end, 
it is convenient to first discretize Eq. ( pT]) 



dt 



= H + J2 M^'ji^) - Ht)) + V^{h). (22) 



where Jij in Fourier space has the form 



+ i.o(p'+'z'), 



(23) 



where A = [LqMI j . The infinite range model is the 
same as in the elastic interface problem 



dt 



H + J{h-h^it))+r],{h), 



(24) 



where h = J^i ^i/N, J = A+Vq and A'' is the system size. 
The mean-field behavior depends on the shape of the ran- 
dom potential: for cusped potentials one obtains that the 
velocity of the interface grows linearly for H > He 



(25) 



A complete mean-field analysis, including the form of 
response and correlation function, can be found in 
Refs. [HJUl. 



To g o beyond mean-field theory, Narayan and Fisher 
|p6| , ^4| have devised a functional renormalization group 
scheme that allows one to obtain the value of the upper 
critical dimension and an estimate of the scaling expo- 
nents. Their method is based on an expansion around 
mean-field theory, using the formalism of Martin-Siggia- 
Rose. They construct a generating functional for the 
response and correlation functions, introducing an auxil- 
iary field h{x,t), 



Z = / {dh){dh)e^p i. i / d^-'x dthF{h,T]) 



(26) 



where 



F{h,ii) = 



dh{x, t) 
dt 



iyoV^h{x,t)- / d'^-'^x' K{x-x'){h{x,t)-h{x,t))-T]{x,h)-H. 



(27) 



Following Ref. |2q], we introduce a new field 



(28) 



which represent of coarse grained version of h, and a cor- 
responding auxiliary field (f). After averaging over the 
disorder one obtains an effective generating functional 



z= {d4>){d4>)eMS{4>,<t>)), 



(29) 



whose saddle point value corresponds with the mean-field 
theory. Narayan and Fisher carried out an expansion 
around the saddle point to obtain the correction to the 
mean-field theory. In our problem everything works like 
in Ref. [M , the only difference being in the form of the 



interaction kernel 
tion is in our case 



q in |g6[). The eflfective ac- 



^d—l J An?, 
£• = I d'^-^x dtH4){x, t) + / ~-uj){-iuj + . " 



d'^^^x dtdt'(j){x, t)C{vt - vt' + (j){x, t) - (f){x, t'))(j){x, t), 



(30) 



where the function C{x) is the mean-field correlation 
function. Other terms resulting from the expansion 
around the saddle point can be seen to be irrelevant. 

To obtain the upper critical dimension, we rescale 
space and time x = bx', t = Vt' , 4> = b'^<l>' , 4> = b^^'^^^cp' 
and H = b~^^'^H', requiring that the Gaussian part of 
the action remains invariant. Simple power counting 
gives 



For d > 3 all non linearities decay to zero at large length- 
scale and the theory is Gaussian, while for d < 3 an infi- 
nite set of non linear terms becomes relevant. The upper 
critical dimension for this problem is therefore dc = 3. 
This result differs from the one obtained for elastic in- 
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terfaces, for which dc = 5j_but agrees with the result 
for contact Hne depinning ||46|] . The similarity between 
the two problems lies in the non-local kernel that scales 
linearly with the momentum at long length scales. 

In order to apply these results to the experiments 
we have to make sure that the linear part of the ker- 
nel dominates on the length scales of interest. Long- 
range effects become relevant for length scales larger than 
L ^ 2i:vq / ^qM"^ . In typical ferromagnets, fig Ms ^ 1 and 
1^0 ~ 10"^ (in IS units) (see page 713 of 0). This im- 
plies L ~ 10^^ 10~® m, which is of the order of the do- 
main wall thickness. From this calculation we conclude 
that the effect of the surface tension can be neglected 
with respect to the long-range kernel. 

Above the upper critical dimension mean-field results 
are valid, while for d = dc we expect logarithmic correc- 
tions. To obtain the value of the exponents below the 
upper critical dimension one should perform a functional 
renormalization group along the lines of Refs. ]25| , ^ , ^ . 
This has been done in Ref . ^] in the case of a kernel scal- 
ing linearly in momentum space. However, in many ex- 
perimental situations the magnetization is perpendicular 
to the plane of the film and this analysis does not apply. 
The issues of Barkhausen effect and domain growth in 
thin films deserve further investigations that are beyond 
the scope of this paper. 



V. CRITICAL EXPONENTS FOR CONSTANT 
APPLIED FIELD 

In this section we derive scaling relations between 
the exponents that characterize the depinning transi- 
tion. When the external field is increased monotoni- 
cally and adiabatically the interface moves in avalanches 
of increasing size. The exponents describing avalanche 
distributions can be compared with experiments on the 
Barkhausen effect. We have to keep in mind that most 
experiments are performed with a non zero applied field 
rate in presence of demagnetizing field. We expect how- 
ever that the distributions dX H = He should scale as in 
the case c ^ and k ^ 0. 

The avalanche size distribution close to the depinning 
transition scales as 



P{s) ^ s~^f{s/so), 



(32) 



where the cutoff scales as sq ^ {H - Hc)-^/" and IS re- 
lated to the correlation length ^ by 



So ~ C 



(33) 



where C is the roughness exponent (Fig. |^) . The correla- 
tion length diverges at the depinning transition as 



e^(i^-i^,)-^ 



(34) 



- = j/(rf-i + C). 

a 



(35) 



The average avalanche size also diverges at the transition 
{s)^{H-H,)-\ (36) 
where 7 is related to r and a by 
(2-r) 



7 = 



(37) 



An additional scaling relation can be obtained consider- 
ing the susceptibility |2^] which is proportional to (s) and 
scales as 



d{h) 

Ilh 



This relation together with Eq. |37j implies 

T = 2- 



1 + K 



v{d-l + 0' 



(38) 



(39) 



The other exponent relevant for the Barkhausen effect 
describes the distribution of avalanche durations 



P[T)^T—g{T/To), 



(40) 



where the cutoff diverges at the transition as Tq ^ 
{H - H^)-^/^. From Eq. (||) and the relation Tq - 
we obtain a = l/zu and 



a = 1 



u{d- 1) - 1 



(41) 



We note that all relations ( p3[ - [4l| ) are valid also for other 
interface problems provided d < dc- 

For our case in d = 3, which corresponds to the upper 
critical dimension, we have C = 0, z = 1 and = 1, that 
inserted in the previous expressions give r = 3/2 and 
a = 2. These exponents agree very well with experimen- 
tal results in the limit of adiabatic driving (c — > 0) . More- 
over, we obtain that the average avalanche size scales 
with the duration as 



{S{T)) - 



(42) 



which implies 



which has been recorded experimentally in Ref. Q. It is 
interesting to compare these results with the exponents 
obtained for three dimensional elastic interfaces. In that 
case the e-expansion gives C = 2/3, z = 14/9 and v = 3/A 
which imply r = 1.25 and a — 1.43 p5| , p6| , ^ . Simula- 
tions give slightly different values, r ~ 1.3 and a ~ 1.5 
pTf . In any case, the values are significantly lower than 
the experimental results. 

When the experiment is performed in absence of de- 
magnetizing fields, as for example in frame geometries, it 
is possible in principle to measure the exponent close to 
the depinning transition. In this regard, several exper- 
iments, discussed in section II, support the mean-field 
prediction v ^ {H — He). Vergne et al. M have observed 
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the growth of the size of the Barkhausen jumps as the 
field is increased. From a measurement of this kind it 
should be possible to obtain an estimate of the exponent 
7. We believe that similar experiments are crucial to 
confirm the presence of a depinning transition. 

Finally, we discuss the properties of the power spec- 
trum of the velocity signal. A similar analysis, in the 
context of fiux line depinning, is reported by Tang et al. 
Q . The height autocorrelation function scales as 

{hir, t)h{f',t')) ^\t- tfi/^f{\r- r'\/\t - /^). 

(43) 




FIG. 5. The interface moves between two pinned configu- 
ration in an avalanche of size s ~ Z''"^"''^ 



The scaling of the velocity autocorrelation function is 
obtained deriving Eq. ( ^ ) with respect to time, which 
gives a power law decay with exponent 2{C,/z — 1). The 
power spectrum of the velocity signal at some fixed space 
location r scales therefore like 

S,{io)^u;^; 7^ = l-2C/z. (44) 



When the velocity is averaged over the whole system we 
expect instead 

S4lo)^J; ^ = l-{2C + l)/z. (45) 



In mean-field theory — 0, which implies = 1 and 
^ = 0. It is interesting to compare these results with 
three dimensional elastic interfaces for which ~ 0.1 
and ijj = —0.6 The direct comparison of these val- 
ues with experimental results is not straightforward due 
to the complexity of the measured spectra. We expect 
the exponents derived from the depinning transition to 
describe the low-frequency part of the power spectrum, 
while for high frequencies we observe a 1//^ decay [ p8[ . 
For low frequencies experiments find exponents ranging 
from 0.5 to 1. This range of value lies between the pre- 
dictions for -0 and ip. We considered the possibility of a 
crossover effect, since in the typical experiment the pick- 
up coil is much smaller than the system size. Depending 
on the domain structure and the coil size the experimen- 
tal exponents could lie anywhere between the averaged 
and non- averaged result. We tested experimentally this 
hypothesis, varying the size of the pick-up coil but we 
noticed no changes in the low frequency part of the spec- 
trum. To obtain a complete explanation of the power 
spectrum, we should probably take into account the pres- 
ence of many interacting domain walls and magnetic after 
effect. 



VI. DRIVING RATE AND DEMAGNETIZING 
FIELD 

In this section we study the effect of the driving rate 
and the demagnetizing field on the dynamics of the 
model. We study here the infinite range model, that in 
d = 3 should have the same critical behavior of the long 
range model, but it is much simpler to analyze. 

As discussed in , the demagnetizing field has the ef- 
fect of keeping the interface close to the depinning transi- 
tion. We will show that the intensity of the demagnetiz- 
ing field is a relevant parameter controlling the avalanche 
characteristic length. Criticality is reached only when 
this parameter is vanishingly small. A finite driving rate 
changes continuously the critical exponents, as in the 
ABBM model We will numerically show that the 

infinite range model reproduces the results of the ABBM 
model and we will present an argument explaining the 
reason for this behavior. This observation explains the 
success of the ABBM model in fitting experimental data. 

The dynamics of the infinite range model is described 
by the following equation 

^^ = H{t)-kh + jCh-h,{t))+rj,{h), (46) 

where the external field H{t) is increased at constant 
rate and the demagnetizing field Hd = —kh has been 
included. 

To show the equivalence with the ABBM model, we 
sum over i both sides of Eq. ( ^6|) and obtain an equation 
for the total magnetization m 

^ = a-km + Y,V^ih), (47) 

2=1 

where the time dependence of the field is now explicit. 
This equation has the same form of the ABBM model 
provided we can interpret r]i as an effective pinning 
W{m), with Brownian correlations. When the interface 
moves between two pinned configuration W changes as 
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W{m') - W{m) = ^ A77i 



(48) 



where the sum is restricted to the n sites that have effec- 
tively moved (i.e. their disorder is changed). The total 
number of such sites scales as n ^ l'^^^ and in mean-field 
theory is proportional to the avalanche size s = \m' — m\ 
(since s ~ and C = 0). Assuming that the Arji are 

uncorrelated and have random signs, we obtain a Brow- 
nian effective pinning field ||50| 



{\Wim') - ly(m)p) = D\m' - m\, 



(49) 



where D quantifies the fiuctuation in W . The Brownian 
pinning field, observed experimentally in SiFe alloys and 
used in the ABBM model to describe the motion of the 
domain wall, is not due to a long-range correlated disor- 
der present in the material. It is instead the result of the 
collective motion of the interface and therefore represents 
only an effective description of the disorder. 

The main predictions of the ABBM model can be ob- 
tained as follows. We derive Eq. ( ^7| ) with respect to time 
and define v = dm / dt 



dv 
'dt 



c — kv + vf{m), 



(50) 



where f(m) = dW/dm is an uncorrelated random field. 
Expressing Eq. Fol as a function of v and m only 



dv 
dm 



k + f{m), 



(51) 



we obtain a Langevin equation for a random walk in a 
confining potential, given by E(v) — kv~ clog(w). In the 
limit of large m, v is given by the Boltzmann distribution 

P{v,m^ oo) - exp{-E{v)/D) = v" cxp{-kv/ D), (52) 

where c = c/D. 

The distribution in the time domain is obtained by a 
simple transformation and it is given by ]l5[ | 

r,/ N N k'^v'^~^exp(~kv/ D) . 
P{v) ^ Piv, t^oo)^ ^.r(^) ■ (53) 

Eq. (|5^) predicts that the domain wall moves at constant 
average velocity (v) — c/k. The relative fluctuations of 
the velocity diverge in the adiabatic limit c — > 



(54) 



This divergence is due to the singularity at low veloc- 
ities of Eq. ( p3| ) and reflects the presence of a depinning 
transition. For c < 1 the velocity distribution is a power 
law with an upper cutoff that diverges as fc — *■ 0. In this 
regime, the domain wall moves in avalanches whose size 
and durations are also distributed as power laws. The 



avalanche size distribution is directly related to the dis- 
tribution of first return times of a random walk in the 
confining potential E{v). Using scaling relations, it has 
been shown |^ that the avalanche exponents scale as a 
function of c as 



T = 3/2-c/2 a = 2-c. 



(55) 



in agreement with experimental results. 

The scaling of the cutoff of the avalanche distributions 
can be obtained as follows. For k = 0, the cutoff in the 
size distribution scales with H as sq {H-Hc)-'^/'' , and 
similarly for the distribution of durations. When fc > 0, 
the interface experience an effective field H — kh which 
keeps it on average below the depinning transition. We 
assume that the distance from the critical point He is of 
the order of 



AH ^ H - Hr 



kAh, 



(56) 



where Ah is the average variation of the height corre- 
sponding to a variation AH in the field. Since Ah ^ 
{s)AH, the average avalanche size scales as 1/fc, which 
implies 



So 



(57) 



Using similar arguments we can also show that the cutoff 
of avalanche durations scales as Tq ~ k~^ in mean-field 
theory. These results do not agree with the experiments 
presented in section II. We will show in the next section 
that they are a peculiarity of mean-field theory and are 
not obeyed by the equation in c? = 3. 

Finally, we note that avalanches are observed only for 
small driving rates (c < 1). For c > 1 the motion is 
smoother with fluctuations that decrease as c increases, 
in agreement with experiments |15[. 



VII. SIMULATIONS 
A. The infinite-range model 

We simulate the inflnite-range model in order to con- 
flrm its equivalence with the ABBM model. We first 
integrate numerically Eq. (^), using the Runge-Kutta 
method and a random potential composed by parabo- 
las with cusp singularities . We study the velocity 
signal as a function of the driving rate c, and find that 
on increasing c, the dynamics crosses over from avalanche 
dominated motion at low c to a smoother motion at larger 
c (see Fig. ||). We are able to integrate the model only for 
relatively small values of N; therefore it is not possible 
to observe the scaling of avalanche distributions, which 
appear to be dominated by finite size effects. 

We then introduce an automaton version of the infi- 
nite range model, which can be simulated for much larger 
system sizes, and study it for different values of c and k. 
From the results of the ABBM model, we expect that 
the velocity distribution is described by Eq. (m). In the 
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limit c — !■ 0, the cutoff in the exponential is vq = k/D. 
We extract vq from the velocity distribution (see Fig. 0a) 
and we plot it for different values of k in Fig. As 
expected, we observe a linear decay and we find a value 
1/D = 1.3 ± 0.1. We then compute the avalanche size 
and duration distribution in the c limit as a function 
of k. The data collapse perfectly (see Figs. ||and ^ using 
the scaling forms predicted in the previous section 



P(s,fc)^s-3/V(sA;2), (58) 



P{T,k)^T-^g{Tk). (59) 



Next, we simulate the model as a function of c and 
find scaling exponents that depend linearly on the driving 
rate. The avalanche size distribution shows a power law 
for more than four decades. Therefore, it provides a reli- 
able estimate of c, using the relation t — 3/2 — c/2. We 
compute r from the distribution as function of c and ob- 
serve a linear behavior r = 3/2 — c/ (2D), with \/D ~ 1.2, 
which is consistent with the scaling of the cutoff of the 
velocity distribution (Fig. 0b). The value of c obtained 
above can then used to fit the velocity and avalanche du- 
ration distributions and the results are consistent with 
the theory (see Figs. ^ and pi] ). 

Finally, we compute the power spectrum for different 
values of c. We observe a 1//^ decay at large frequency 
and a constant part at low frequencies. The crossover 
point scales linearly with c as in experiments . 
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FIG. 6. The velocity of the interface as a function of time, 
for different values of c. The data have been obtained in- 
tegrating the equation of motion with cusped potential for 
N = 256. 
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FIG. 7. (a) The distribution of velocities in the infi- 
nite-range automaton model as a function of k for A*' — 32696, 
c = 0. (b) The scaling of the 1/vo cutoff with k. The line is 
a fit with slope 1/D ~ 1.3. 
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FIG. 8. (a) The distribution of avalanche sizes in the infi- FIG. 10. The distribution of velocities in the infinite-range 

nite-range automaton model as a function of k for — 32696, automaton model as a function of c for A'^ = 32696, 

c = 0. A line corresponding to t = 3/2 is plotted for com- k = 0.0075. The lines are the theoretical predictions 1 — c. 
parison. (b) The corresponding plot, using scaled variables, 
showing excellent data collapse. 
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FIG. 9. (a) The distribution of avalanche durations in 
the infinite-range automaton model as a function of k for 
N = 32696, c = 0. A line corresponding to a = 2 is plot- 
ted for comparison, (b) The corresponding plot, using scaled 
variables, showing excellent data collapse. 
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FIG. 11. (a)The distribution of avalanche sizes in the in- 
finite-range automaton model for different driving rates for 
N = 32696, k = 0.0075. The fit of the power-law part yields 
T = 3/2- c/2, with c = c/D and 1/D ~ 1.2. (b) The corre- 
sponding distribution of avalanche durations. The power-law 
part is fit with an exponent a = 2 — c. 



B. The long-range model 

A numerical integration of Eq. ( ^ ) poses serious nu- 
merical problems due to the presence of long range non- 
local kernel. Therefore, we study an automaton version of 
the model, which should belong to the same universality 
class. In the automaton model the height is discretized 
and the local velocity can assume only the values v = 0,1. 
For each configuration of the system, we compute the lo- 
cal force according to Eq. (|2l]). Periodic boundary con- 
ditions are imposed on the lattice and therefore we must 
sum the non-local kernel over the images as discussed in 
Ref. To model the disorder, we associate to each 

site on the interface a random number chosen from a 
Gaussian distribution. 

When the local force on a site is larger than zero, 
the corresponding height is increased by one unit and 
we chose a new value for the disorder. Care must be 
taken in choosing the values of the parameters, in order 
to avoid instabilities present in the discretization of the 
kernel ||. 

We consider first the case fc = to confirm the pre- 
dictions about the upper critical dimension. We increase 
the external filed adiabatically up to He (i.e. when the 
interface is pinned we increase the external field until the 
most unstable site reaches the threshold for movement), 
and we compute the integrated avalanche size distribu- 
tion. This distribution scales as 

P^nt{s)= / dHPis,H)^ s-^^+''\ (60) 
Jo 

which yields a s^^ decay in mean-field theory. Similarly, 
for the integrated duration distribution we find a 
decay. The simulation results confirm the predictions of 
the theory (see Fig. p^ . 

Next, we study the model in the adiabatic limit (c 0) 
as a function of k. We compute the distribution of veloci- 
ties (Fig. |l^) and avalanches sizes (Fig. ^) and durations 
(Fig. |l5| ) as a function of the demagnetizing field k. The 
scaling exponents are in agreement with the results of 
the depinning transition in the mean- field, r = 3/2 and 
a = 2. 

However, the scaling of the cutoff of the distributions 
does not agree with the predictions of the ABBM model. 
We find instead sq ^ k~^ and Tq ^ This behavior 

persists in simulations performed at c > 0, where the ex- 
ponent T and a still scale with c as in the ABBM model. 
To obtain a good data collapse, the scaling functions in 
Eqs. ( |58| ) and (|5^) have to be replaced by 



Pis,k)^s-^^^gi{sk), (61) 

P(r,fc)^T-2g2(Tfci/2), (62) 

which are the scaling forms obtained experimentally (see 
section II). 

The precise reason for these results is still not com- 
pletely clear. Recent simulations of a model similar to 
ours, studied in the context of dry friction, suggest that 
the effective pinning field for the long-range model is not 
Brownian |Q . In Ref. ||5^ the cut-off of the distributions 
was related to the shape of the force distribution, but it 
is not clear if this analysis can be applied directly to our 
case, due to the different driving mechanism employed in 
Ref. H. 
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FIG. 12. The integrated avalanche size distribution in the 
long-range automaton model for fc = and L = 61. A line 
with slope —2 is plotted for reference. 
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FIG. 13. The velocity distribution in the long-range au- 
tomaton model for c = and L = 61. A line with slope —1 

is reported for reference. In the inset we show the linear-log 
plot of the same distribution in order to show the exponential 
cutoff. 
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FIG. 15. a) The avalanche duration distribution as a func- 
tion of k in the long-range automaton model for c = and 
L = 61. A line with slope —2 is reported for reference, b) The 
corresponding plot, using scaled variables, showing excellent 
data collapse. 



FIG. 14. a) The avalanche size distribution for c = as 
a function of k for the long-range automaton model with 
L = 61. A line with slope —3/2 is reported for reference, 
b) The corresponding plot, using scaled variables, showing 
excellent data collapse. 



VIII. DISCUSSION 

In this paper we have studied the dynamics of a flexible 
domain wall as it moves through a disordered medium. 
We have derived an equation of motion taking into ac- 
count the effect of different energetic contributions. A 
crucial role is played by dipolar interactions that give 
rise to a demagnetizing field and to a long-range inter- 
action kernel. In absence of demagnetizing field, the do- 
main wall shows a depinniiig transition as a function of 
the field. The long-range interaction kernel set the upper 
critical dimension to = 3, so that mean-field scaling 
should describe the experiments on the Barkhausen ef- 
fect. 

The predictions of the present theory compare well 
with the distribution of Barkhausen jump durations and 
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sizes and with the velocity distribution. In particular, 
we discuss the Unear dependence of the exponents on the 
field driving rate P,p[ |3l| , |3^ and the scaling of the cutoff 
with the demagnetizing field. The agreement between 
theory and experiments is in both cases quantitative. In 
toroidal geometries, when the demagnetizing field is zero, 
we predict a linear dependence of the domain wall veloc- 
ity on the applied field, in agreement with several exper- 
iments on soft ferromagnetic materials pTf . 

We show that the phenomenological model introduced 
by ABBM |l^ is equivalent to the infinite-range domain 
wall. The Brownian correlated random pinning field used 
in [ p^ and experimentally observed in Q is shown to 
arise in the effective description of the motion of center 
of mass of the domain wall. This result clarifies the origin 
of the correlated disorder which could not be explained 
as a simple result of the correlations between the impu- 
rities . While the infinite range model, and therefore 
the ABBM model, quantitatively explains many features 
of the Barkhausen effect, it does not give the correct de- 
pendence on the demagnetizing field, which is instead 
provided by the complete three-dimensional description. 

The power spectrum of the Barkhausen noise does not 
show a marked universality and therefore can not be com- 
pletely explained by our approach. In particular, we ob- 
tain a 1//^ decay at large frequencies, which has only 
been observed in experiments with a single domain wall 
1^ . Other experimental results seem to suggest that the 
exponent changes when the number of domain walls in- 
creases 1^,^,^. Moreover, magnetic after-effect [|ll[ and 
fltix propagation could also affect the results. To obtain 
a quantitative explanations of these results, one should 
analyze the dynamics of many coupled domain walls. 

The presence of many domain walls should affect the 
power spectrum, but not the avalanche distributions. 
When a domain wall starts to move, the demagnetiz- 
ing field increases, creating a larger pinning force on the 
other walls. Therefore, on short time scales the interac- 
tions between the walls are irrelevant. For this reason, 
the avalanche distribution for a single domain wall agrees 
with experiments performed with many domain walls. 

With our approach we can address several other is- 
sues raised in the literature about the Barkhausen effect. 
The partial reproducibility of the Barkhausen signal ob- 
served in recent experiments [|l^,^,^ is explained by 
the quenched nature of the disorder. Pushing the wall 
back and forth through the same disordered region of the 
sample results in the same signal. Deviations from this 
ideal behavior can be expected due to small variations on 
the initial conditions, thermal effects or differences in the 
driving rate. To understand these features it is crucial 
to consider a flexible domain wall instead of a rigid wall 
| p5[ , for which always perfect reproducibility is expected. 

The recent theoretical revival in the study of the 
Barkhausen effect is mostly due to the claim of Ref. Q 
that this phenomenon is an example of self-organized 
criticality (SOC) This claim was challanged in 

Ref. which, based on the results obtained for the 



RFIM, concluded that scaling in Barkhausen effect is due 
to the presence of a "plain old" critical point. The ques- 
tion concerns the origin of the cutoff in the power-law 
distributions. According to the analysis of Refs. [^-^, 
the cutoff would be determined by the variance of the 
random-fleld distribution. As far as we know, no experi- 
mental evidence of a critical point of this kind has been 
reported in the literature. 

We have experimentally observed that the cutoff of the 
distributions is determined by the demagnetizing field, in 
agreement with our theoretical analysis. In our model, 
the critical point is reached only by fine tuning to zero 
the driving rate and the demagnetizing field, performing 
the limits c — > and fc — > in the given order. This 
result implies that the idea of a critical point without 
tuning parameters does not apply to the Barkhausen ef- 
fect. It is interesting to remark, however, that the picture 
revealed by our approach is similar to the one presented 
in Ref. where it was shown that criticality in SOC 
models arise by the fine tuning to zero of the driving 
rate and other parameters. In light of these results, one 
should review the definition of SOC ||5^, dropping the 
notion of "criticality without fine-tuning" . In this re- 
stricted sense, the Barkhausen effect can be considered 
as a good example of SOC. 

The present approach to the Barkhausen effect, based 
on the depinning of a ferromagnetic domain wall, ap- 
plies to three dimensional soft ferromagnetic materials, 
which arc frequently used in experimental studies of the 
Barkhausen effect. For hard ferromagnet and rare earth 
materials, where strong local anisotropies prevent the for- 
mation of straight domain walls, a different approach is 
needed. Disordered spin models like those presented in 
Refs. |20|-p^ seem more appropriate. We did not discuss 
here the issue of domain nucleation and growth in thin 
films (two dimensional ferromagnets). Depending on the 
material properties and the sample geometry, the domain 
walls are either fractal or self-affine as in our case. In the 
second case, we expect that the framework of the depin- 
ning transition could be relevant. 
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